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ABSTRACT 

A number of recent lines of evidence point towards the presence of hot, out- 
flowing plasma from the central regions of compact Galactic and extragalactic 
X-ray sources. Additionally, it has long been noted that many of these sources 
exhibit an "excess" continuum component, above ~ 10 keV, usually attributed to 
Compton Reflection from a static medium. Motivated by these facts, as well as by 
recent observational constraints on the Compton reflection models - specifically 
apparently discrepant variability timescales for line and continuum components 
in some cases - we consider possible effects of out-flowing plasma on the high- 
energy continuum spectra of accretion powered compact objects. We present a 
general formulation for photon downscattering diffusion which includes recoil and 
Comptonization effects due to divergence of the flow. We then develop an ana- 
lytical theory for the spectral formation in such systems that allows us to derive 
formulae for the emergent spectrum. Finally we perform the analytical model 
fitting on several Galactic X-ray binaries. Objects which have been modeled with 
high-covering-fraction Compton reflectors, such as GS1353-64 are included in our 
analysis. In addition, Cyg X-3, is which is widely believed to be characterized 
by dense circumstellar winds with temperature of order 10 6 K, provides an inter- 
esting test case. Data from INTEGRAL and RXTE covering the ~ 3 — 300 keV 
range are used in our analysis. We further consider the possibility that the widely 
noted distortion of the power-law continuum above 10 keV may in some cases be 
explained by these spectral softening effects. 
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Subject headings: accretion — stars: radiation mechanisms: nonthermal - stars: 
black holes - stars: individual (GS1353-64, Cyg X-3, Cyg X-l, GX339-4 ) 

1. INTRODUCTION 

Recent observational and theoretical evidence suggest that accretion-powered X-ray 
sources, of both the Galactic and extragalactic variety, may exhibit outflowing plasma, i.e. 
winds, emanating from a compact region near the central source (e.g. Elvis 2003; Arav 2003; 
Brandt & Schulz 2000; Proga & Kallman 2002). Comptonization effects in those putative 
outflows are likely to alter the intrinsic continuum spectra of accretion powered compact 
objects. The basic idea is that electron scattering of photons from a central source entering 
the expanding outflow experience a decrease in energy (downscattering). The magnitude of 
this decrease is of first order in v/c and in E/m e c 2 where v is the outflow speed, c is the 
speed of light, E is the initial photon energy, and m e is the electron rest mass. 

The basic idea is depicted in Figure 1. There present a simple explanation of the 
diverging flow effect on the photon propagation through the medium. A photon emitted 
outwards near inner boundary and then scattered at a certain point by an electron moving 
with velocity v 1; is received by an electron moving with velocity v 2 as shown with frequency 
v 2 = V\ [1 + (v x — v 2 ) • n/c] where n is a unit vector along the path of the photon at the 
scattering point. In a diverging flow (vi — v 2 ) • n/c < and photons are successively red- 
shifted, until scattered to an observer at infinity. The color of photon path (in Figure 1) 
indicates the frequency shift in the rest frame of the receiver (electron or the Earth ob- 
server). On the other hand, referring to the right-hand side of Figure 1, in a converging flow 
(v x — v 2 ) • n/c > and photons are blueshifted. 

The classical Compton (recoil) effect < AE > /E ~ —E/m e c 2 has been well understood 
for a long time. Basko, Sunyaev & Titarchuk (1974), hereafter BST74 first studied the 
downscattering effects in the interaction of X-ray radiation of the central source with the 
relatively cold atmosphere of the optical companion in a binary system. They predicted 
the shape of the X-ray reflection spectrum of the companion and applied these results to 
the Her X-l system. X-ray observations of Her X-l by a number of groups (e.g. Sheffer et 
al. 1992, Still et al. 2000 and etc.) confirmed this prediction. In particular, the numerical 
calculations by BST74 demonstrated that the reflection of the bremsstrahlung spectrum 
distorted the high-energy continuum above 10 keV producing the characteristic continuum 
excess feature, or "bump", in the > 10 keV spectrum. This feature was also identified by 
Sunyaev & Titarchuk (1980), hereafter ST80, who found that the transmitted, downscattered 
spectrum is formed as a result of the reprocessing of X-ray radiation from the central source 
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in an ambient spherical cloud. From these facts one can conclude that this distortion of the 
continuum is not an intrinsic property of the particular space photon distribution, but it is 
rather a result of effective downscattering (in terms of the mean number of scattering and 
the energy loss per scattering) X-ray photons diffusing through the intervening medium. 

In this paper we show that this is precisely the case. We formulate the problem of the 
photon diffusion in generic terms and demonstrate a specific solution for bulk motion where 
effects of a recoil and divergence of the flow are taken into account. The emergent spectrum 
as an outcome of this solution has all these aforementioned features of the downscattering. 

The details of the radiative transfer problem taking into account downscattering effects 
(recoil and Doppler effect in the divergent flow) and its solution are given in §2. In section 
3, we apply our model to observational data for several Galactic X-ray binaries. We attempt 
to demonstrate the Compton downscattering effects in relatively cold outflows for which 
temperature is of order 10 6 K and disk atmospheres which are necessary constituents of X- 
ray sources. For example, Cygnus X-3 is believed to be characterized by dense circum-stellar 
(or circum-disk) winds, provides an interesting test case. A number of additional sources, 
which have been noted in the literature to exhibit evidence for strong "Compton reflection", 
e.g. GS1353-64 and GX339-4, have also been included in our analysis. We consider the 
possibility that this well documented distortion of the power-law continuum above 10 keV, 
may in some cases be due to downscattering from out-flowing plasma rather than from a 
static reflecting media. We further speculate that this could possibly apply to extra-galactic 
(i.e. AGN) as well as Galactic accretion objects (§4), although we do not present any analysis 
of AGN here. A summary and conclusions follow in §5. 



2. RADIATIVE TRANSFER IN A BULK OUTFLOW 

2.1. The Basic Downscattering Problem 

The problem of photon propagation in a fluid in bulk motion has been studied in detail in 
a number of papers [see e.g. Blandford & Payne (1981), hereafter BP81; Payne & Blandford 
(1981), hereafter PB81; Nobili, Turolla & Zampieri (1993); Titarchuk, Mastichiadis & Kylafis 
(1997); Titarchuk, Kazanas & Becker (2003, hereafter TKB03); Laurent & Titarchuk (1999), 
(2001), (2004, hereafter LT04]. In particular, TKB03 present a general formulation and a 
solution of the spectral formation in the diverging outflow. They demonstrated that the 
resulting spectrum can be formed as a convolution of energy and space diffusion solutions. 
The spread function, given by a Green's function formulation, applied to a monochromatic 
injection yields a redshift-skewed line which is a power law at low energies. Furthermore, 
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Monte-Carlo simulations show that the line width has a strong dependence on v/c, optical 
depth r and the energy of the monochromatic line E only when the flow temperature kT e < 1 
keV (LT04). TKB03 and LT04 also establish that the shape of the spectral line should be 
closely related to the photon source distribution in the flow. Comparison of the theoretical 
spectra with data show that observed red-skewed lines are formed as a result of transmission 
of the X-ray radiation through outflows of moderate Thomson optical depth (tqZ 1)- 

To apply these Radiative Transfer results to observational data, we have developed a 
generic analytical formulation which leads to a simple analytic expression for modification 
of the emergent spectrum due to recoil and velocity divergence effects in the flow. For the 
recoil effect, we extend the results of ST80 who show that the downscattering feature, or 
"bump", can appear superposed on power- law spectra with energy spectral indices a < 1 as 
a result of photon diffusion through a static cloud. From ST80 and calculations presented 
here, one might conclude such bumps are not a necessarily a feature of disk reflection, but 
could instead be a generic feature of the photon reprocessing in a relatively cool, ambient 
plasma characterized by temperatures of order 10 6 K. 



Let N(r) = N (r /r) b be the radial number density profile of an outflow and let its 
radial outward speed be 



is the proton mass). The Thomson optical depth of the flow from some radius r to infinity 
is given by 



where N e (r) = N(r) is the electron density, <tt is the Thomson cross section, r is a radius at 
the base of the outflow. Because our final results are independent of the velocity and density 
profiles below for simplicity of presentation we use b = 2 (a constant velocity outflow). In 
this case r = r(r ) = (T T A^or . 

The transfer of radiation within the flow in space and energy is governed by the photon 
kinetic equation (BP81, Eq. 18) for the photon occupation number n(r, u), which in steady 



2.2. Main Equations and Solution of Downscattering Problem 





(2) 
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state reads 



v b/c _ 1 _ , 1_ N 1 V • (vb/c) On 1 d(z 4 n) .. . , oX 
Vn+— V-(-Vn) + - L^L^ + v ; =-j(r,z/), (3) 



/c 3k k 3 k dz z 2 dz 

where z = hu/m e c 2 is the dimensionless photon energy, k = N e {j-)a^ is the inverse of the 
scattering mean free path, Vb = vi,e r , is the flow velocity, e r is the radial unit vector and 
j(r, v) is the photon source term. 

The spectral flux F(r, v) (PB81) is given in terms of n(r, v) by 

and must satisfy the following boundary conditions: (a) Conservation of the frequency inte- 
grated flux over the outer boundary, namely 

L^n = J F(r, v)dv oc r~ 2 as r — > oo. (5) 

(b) In terms of the specific intensity of the radiation, the second inner boundary condition 
specifies that the escaping radiation from the inner boundary be equal to the radiation 
entering the wind shell through the inner boundary (see Fig. 1 for the wind geometry). 
We neglect here the possible absorption of the radiation at the central source (black hole). 
Thus the occupation number should satisfy some kind of the reflection condition at the 
inner boundary r = r . Namely the photon flux at the inner boundary is zero (all photons 
scattered in the outflow shell and escape through the inner boundary subsequently return). 

L^n(r ) = 0. (6) 

In this formulation of the problem the number of photons emitted in the wind shell equals 
to that escaped to the Earth observer (i.e. the photon number is conserved). This can be 
proven using Eqs. (3), (4-5). Consequently, if the high energy photons lose their energy in 
the way out but the number of photons is conserved it has to be the accumulation of the 
photons at a particular lower energy band. This accumulation effect was previously noted 
by ST80 for the case when the photon energy loss (downscattering) was determined by the 
recoil effect only (see Fig. 10 in there). In section 2.3 we demonstrate this accumulation 
effect when the recoil and flow divergence dowscattering effects are taken into account. 

It is worth noting that in TKB03 the goal was to demonstrate the pure effect of the 
flow divergence on a spectral line as a result of multiple scatterings with the plasma free 
electrons. They thus omitted the recoil term in the left hand side of equation (3). They 
were also interested to applying their solution to the spectral formation of iron lines for 
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which photon energies are about 6.4 keV and lower. For such energies the recoil effect can 
be safely neglected (see ST80). On the other hand if one is interested in the modification 
of continuum spectrum at energies higher than 10 keV by multiple downscattering events 
in the "cold" medium this effect must be taken into account. Below we show how one can 
generalize the TKB03 solution including the recoil term in the radiative transfer equation. 
We apply a particular method for the separation of variables suggested in Titarchuk (1994), 
hereafter T94, and successfully used in TKB03, to obtain our solution. 

In order to further proceed with the solution derivation we should note that Laurent 
and Titarchuk (2004) find that the coefficient in the diverging term of equation (3) e = 
V • (vb/c)//c can be replaced by a constant e. Strictly speaking e is a function of r and 
consequently of r, namely e = 2(i> 6 /c)/r. 

Using the Fokker-Planck equation [see Eq. (3) in TKB03] and a method developed by 
Titarchuk, Mastichiadis & Kylafis (1997) (see Appendix D there) one can relate the mean 
energy change per scattering, < AE >, for photons undergoing numerous scatterings in the 
flow to P — v/c: 

< AE >^ -\fVw/{cK)]E » -(2ff3/r )E (7) 

where v = ve r is the flow velocity, e r is the radial unit vector, k = N e a T is the inverse of 
scattering mean-free-path /. The numerical factor /, in formula (7) is of order unity and 
LT04 obtain its precise value using their MC simulations. 

In LT04 (specifically, Figure 2 of that paper) the escape photon distribution ip(t) for 5 
energy bands is present. In the plot the time is given in light crossing time units t cross = Ar/c 
where Ar is the outflow cloud thickness. Photons which escape without any scattering are 
at 6.6 keV. The model parameters kT e = 0.1 keV, r = 4, (3 = v/c — 0.1. 

The simulated time distribution is fitted by an exponential ip(t) = exp(— at/t cross ), 
with a equal to 0.67 identical to that predicted by diffusion theory [see Sunyaev & Titarchuk 
in 1985 (hereafter ST85)]. ST85 show that the average number of scatterings N av in the 
shell of optical depth r is 3tq/8 so the average photon scattering time is t av = N av l/c = 
3T t cross /8, where / = l/N cr T = Ar/r . ST85 also show that the time distribution for 
scattered photons in any bounded medium is an exponential, tp(t) = exp(— t/t av ) where 
C N = l/t av (J °° f(t)dt = l). 

In this particular case f(t) = Cn exp(—t/t av ) = Cn exp[— (2/3)/(r /4)(t/t cross )], pre- 
cisely what is obtained in the LT04 simulations. Thus, the average energy of photons escap- 
ing after N av — scatterings is < E > sc ~ (1+ < AE >) Nav . Using formula (7) one can obtain 
that 

< E > sc ^ (1 - 2ff3/r ) N ™E . (8) 
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For the particular case of t = 4 and (3 — 0.1 LT04 find that the original photon energy 
E = 6.6 keV is reduced to E Q = 5.4 keV, after N av — scatterings in the outflow [the emergent 
spectrum is shown in Fig 1, (LT04) left panel]. The analytical estimate obtained from 
formula (8) is close to this value of < E > sc = 5.4 keV for / fa 2/3. Thus one considers an 
approximation in which the diverging term in equation (3) is independent of space variable. 



2.3. Downscattering solution: Emergent spectrum 

According to a theorem (T94, appendix A), the solution of any equation whose LHS 
operator acting on the unknown function n(r, is) is the sum of two operators L r and L u , which 
depend correspondingly only on space and energy and the RHS, j(r,v), can be factorized, 
i.e. 

L r n + L u n = -j(r, v) = -f{r)tp(v). (9) 
The boundary conditions independent of the energy u, i.e. 

L^n = as r -> oo, L^n = for r = r , (10) 

n is given by the convolution of the solutions of the time-dependent problem of each operator, 
namely 

poo 

n(r,u) — / P{r,u)X{y,u)du. (11) 



Above u is the dimensionless time (which can be the Thomson dimesionless time «t — N a T ct 
depending on the specific forms of operators L r and L v ) and P(r,u) is the solution of the 
initial value problem of the spatial operator L r 

dP 

~^ = L r P, P(r,0) = f(r) (12) 

with boundary conditions 

4 1} P = as r^oo, L^P = at r = r , (13) 

and X(v, u) the solution of the initial value problem of the energy operator L v 

8X 

— = L V X, X(z,0) = V (z) (14) 

with boundary conditions 

z 3 X -> when z -> 0, oo. (15) 

Thus we have 

Lr P = _Z^. V P+^V-(ivP). (16) 

Hi O Hi Hi 
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Following, the LT04 arguments we replace the coefficient of the diverging term by a 
constant e = 2g(t> b /c)/r T <C 1 where a numerical factor q is order of unity (see details in 
the end of §2.2). 

In the case where the photon energy is due to a recoil effect and Comptonization effects 
in the diverging flow, we have 

f = LuX = i £ + (17) 

Ou 3 Oz z l oz 

X(z,0)=p(z)/z 3 (18) 

with boundary conditions 

z 3 X -> when z -> 0, oo. (19) 

We transform Eq.(17) introducing a new unknown function F = e^ 4£ ^ 3 - )?t z 4 X for which equa- 
tion has a form 

The problem for equation (20) with appropriate initial condition Y(z, 0) = zip(z) and 
boundary conditions Y — > when 2^0, oo is an initial value problem for the first order 
partial differential equation and it can be found using the method of characteristics (see e.g. 
TKB03). The differential equation for the characteristics is 

dz 

du = 77 (21) 

(e/3)z + z 2 V ; 

that solution is 

u = -[\n{z + e/3)/z - \n(z + e/3)/z )\, (22) 
e 

where Zq is the dimensionless energy at u — 0. Because Y(z , 0) = Zoip(zo) is conserved along 
the characteristics the solution of the problem (17-19) is 

J(z,u) = z 3 X{z,u) = e' i£u/3 ^ £ \z,u)]/[z^ £ \z,u)}, (23) 

where z = tp^(z,u) = (e/3)/[(l + e/3z) exp(— eu/3) — 1] is found from Eq. (22). Substitu- 
tion of J(z,u) from Eq.(23) into Eq.(ll) gives us the emergent spectral shape 

1 ru m ax,e{z) 

T E {z,e) = [r- 2 F(r,z)}\ T ^ oc - / e^ 3 [^ e \z,u)]-M^ e) M]V(u)du (24) 

z Jo 

where u maXj£ (z) = (3/e)ln(l + e/3z) 3> 1 and V(u) oc [t~ 2 <9P/<9t[t, u]\ t ^q using the ex- 
pressions for F(t,z) and P(t,u) [see Eq. (4), Eqs. (12-13) for the definition of F(t,z) and 
P(t,u) respectively]. 
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This result is a generalization of ST80 and TKB03 results. ST80 derived the spectra 
when the downscattering effects due to the recoil was taken into account. On the other hand 
TKB03 considered the downscattered spectra as a result of the divergence of the flow. For 
example, one can obtain the ST80 formula (36) for the recoil spectrum assuming e — > in 
Eq. (24). In fact, 

\imip( £ \z, u) = (1/z — u)~ l 

s— »0 

and u maxfi (z) = 1/z. 

That formula (24) fore = (the recoil case) is generic and it is valid for any geometric 
configuration of the plasma cloud, e.g. a disk, or a spherical cloud, as well as for any photon 
source distribution within the cloud, e.g. uniform, or central illumination distributions. In 
fact, what we use here to derive this formula are the particular properties of the diffusion 
operator of the left hand side of the main equation (9), the boundary conditions (10) and 
the source term in right hand side of equation (9). Namely i. the diffusion operator should 
be a sum of two operators, one is the space diffusion operator and another one is the energy 
operator, ii. the boundary condition are independent of energy and iii. the source function 
is factorized [or presented as a linear superposition of the products of f(r)ip(y)]. All these 
mathematical properties are generic for the diffusion problem and independent of any specific 
geometry and spectral and space source distribution (see ST80 and T94 for more details and 
particular examples). 

Below we demonstrate how formula (24) can be simplified by exploiting the fact that 
downscattering energy change AE/E is proportional to e and E/m e c 2 . 

It should be noted that the downscattering modification of the spectrum occurs when 
photons undergo multiple scatterings u. It is always on the order of the average number 
of scatterings N av = J °° uV[u)du/ J °° V[u)du. Thus we can expand the integrand function 

W&(z,u) = e- 4eu/3 [^ (e) ('2 ) M )] _ V[^ (e) (^ ! w)] over u in formula (24): 

W {e) (z, u) « W {£) (z, 0) + Wi £) '(z, 0)u. (25) 

This expansion is valid because we consider the case when N av <^ l/e (tt,o~1, £ 1). 
Substitution of Eq (25) into Eq. (24) gives us 

T[v) oc z~ 1 [W^(z,0) + N av Wjf^(z,0)], (26) 

where z~ 1 W ( ' £ \z, 0) is the incident spectrum (p(z) and 

W^\zM/z = V(z) ({ln[^o)]K k=o^k=o -4s/3j , (27) 
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where z = ip( £ \z, u) = (e/3)/[(l + e/3z) exp(— eu/3) — 1] and 

= *"<! + £). (28) 

Substitution of equations (27, 28) and </?(z) = z^W^^z, 0) into (26) leads us to the formula: 

T E {z,e) oc <p(z){l - (Ae/3)N av + N av z 2 (l + e/3z)[\n(z oV (z ))]'J u=0 }. (29) 

The most interesting case is the downscattering modification of Comptonization spectrum 
(see e.g. ST80 and Titarchuk 1994) that can be well fitted by a powerlaw spectrum with an 
exponential cutoff, namely 

Vcomp(z) ~ z~ a exp(-z/z*). (30) 

The cutoff energy E* is related to the Compton cloud electron temperature kT e , namely 
E* ~ 2kT e . For ip C om P (z), it is evident that 

[ln(z (p(z ))]' Zo \u=o = {l-a)/z- 1/z* (31) 

Using equation (31) we transform formula (29) as follows 

F E (z, e) oc ¥ C om P (z){l - (Ae/3)N av + N av z(l + e/3z)[{l -a)- z/z*]}. (32) 

The modification of the absolute normalization of the spectrum due to the downscattering 
effects is quite obvious — the relative change of the normalization is 1 — (1 + a)eN av /3 . 
Thus we can rewrite Eq. (32) as the final formula of the spectral shape as follows 

T E {z,e) oc ip CO mp(z){l + N av z[(l - a) - e/3z* - z/z*\, (33) 

where N av = N av /[1 — (1 + a)eN av /3]. The second term in parenthesis of formula (33) 
describes the pile up and softening of the Comptonization spectrum (p CO mp(z) due to the 
downscattering effect in the outflow. In Figure 2a we present ETe diagram for various N av 
and e. The downscattering (accumulation) bump and softening of spectrum at high energies 
are clearly seen in this plot [see also Fig. 2 where a ratio of the Comptonization models to 
the incident spectrum (30) is plotted]. 

Below we apply formula (33) to the X-ray spectral data for a number of sources. The 
self-consistency of application of this formula for data fitting can be checked by comparison 
of the best-fit parameter N av and u max ^[z max ) [the formula for u maX;£ (z) is shown just after 
Eq. (24)]. Our inferred best-fit parameter N av u max ^(z max ) for all fits and thus our 
formula (33) as the analytical approximation of (24) is valid for all cases considered. 
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3. DATA ANALYSIS 

To test these ideas we obtained public data from the INTEGRAL and HEASARC 
(RXTE) archives for a sample of Galactic X-ray binaries (Table 1). The Cygnus region 
was covered extensively by INTEGRAL in the early stages of that mission, Performance 
Verification Phase, (November-December 2002), subsets of which are in the public domain 
at the time of this submission. Specifically, we selected coverage of Cyg X-3 and Cyg X-l 
obtained during revolutions 22-25. Statistics are limited by the small amount of useful data, 
given the early problems with data gaps. We extracted spectra from JEM-X and SPI, using 
the OSA release version 3 software. Subsequent model fitting was done with the XSPEC 
spectral analysis package (version 12.10a:) modified to include a local implementation of the 
model described in section 2. We also analyzed the JEM-X and SPI data for Cyg X-l, 
applying the same model. For comparison, we also obtained and analyzed contemporaneous 
RXTE data for both sources. In all cases, we used the RXTE calibration database and 
software release current as of April, 2004. We also utilized RXTE data for the additional 
sources GX 339-4 and GS 1354-63, each of which have been modeled by various other groups 
- e.g. Gilfanov et al (1999); Nowak, Wilms & Dove (2002); Pottschmidt et al (2003); Vilhu 
et al (2003) - as a power-law- exponential plus a Compton "reflection" continuum (see e.g. 
Magdziarz & Zdziarski 1995). For example, GS 1353-64 was found to require a large reflection 
continuum, R ~ 0.3 — 0.6 (Gilfanov et al 2003). In addition to the Comptonized continuum, 
iron line structure - a gaussian line component and absorption edge feature - was found 
to improve the quality of fit at low energies. In addition, a systematic error component 
of 1% was added to the PCA data. We note that there is a significant cross-calibration 
discrepancy between the INTEGRAL SPI and JEM-X instruments (e.g. Paizis et al 2003). 
We have assumed that the absolute flux calibration of SPI is reliable (Attie' et al. 2003; 
Sturner et al. 2003), and renormalized the JEM-X model fits accordingly. 

In Table 1 we list our source sample, and the inferred parameters of the out-flow- 
Comptonization from our analysis. Here r = a + 1 is the (photon) spectral index. N av is 
the average number of scatterings experienced by a typical photon in the outflow (section 
2). The origin of the data RXTE (PCA plus HEXTE) or INTEGRAL (JEM-X plus SPI) 
are also indicated. We note that two Cyg X-3 observations represent very different intensity 
states for that source. For that matter, the GX 339-4 and GS 1353-64 observations used 
correspond to high-intensity states for each of those sources (but in case, the spectral energy 
distributions represent the low-hard state). 

Figure 3 further illustrates the downscattering effect. Plotted there is our best fit model 
curve, using GX339-4 as a test case, compared to a simple exponentially- folded power-law 
fitted to the same data, with an accompanying residual plot. This illustrates the improvement 
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in the fit above 10 keV. It is clear from these plot that the downscattering effects lead to an 
improved fit above ~ 10 keV. Also one can see the photoabsorption features below 10 keV 
and higher K-edge (in 7-9 keV energy band) in the residual plot (Fig. 3a). It is worth the 
similar edge features along with the strong K a lines are detected during X-ray superbursts 
(see Figs. 5 and 9 in Strohmayer & Brown 2002). There is a high probability that they are 
originated in the radiation driven outflows during burst events. 

In Figure 4, we show an example of one of our model fits, in this case to GS 1353-64 
(Figure 5 is the same result, but plotted in photon space). 

4. DISCUSSION 

We have shown that downscattering modification of the primary photon spectrum by an 
outflowing plasma is a possible mechanism for producing the continuum excess in the ~ 10 
keV spectral region. This is usually attributed to Comptonization by a static reflector, such 
as a downward- or obliquely-illuminated accretion disk, although the overall continuum form 
differs from that of the basic Compton reflection form. We thus suggest, that in at least 
some cases, the outflow downscattering effect rather than the standard Compton reflection 
mechanism is responsible for the observed "excess" hard- X-ray continuum. 

It is reasonable to expect that outflowing plasma, in the form of stellar or putative 
disk winds will effect the emergent spectra of compact binaries. Collimated outflows are 
well known in certain objects, the so called Galactic micro-quasars, and there is strong un- 
collimated or weakly-collimated outflow in additional objects; as noted Cyg X-3, where in 
fact spatially extended emission has been resolved (Heindl et al. 2003), and of course the 
stellar winds in Cyg X-l have been extensively studied. Disk winds may have been observed 
directly in Circinus X-l, e.g. Brandt & Schulz (2000). Recently, evidence for disk winds 
in AGN has emerged (e.g. Elvis 2001; Arav 2003). Further observational confirmation is 
needed, but if present, these winds could produce the downscattering effects we propose. 

In addition to evidence for outflows, recent observations suggest that the putative reflec- 
tor is, in at least certain objects, necessarily much larger than any reasonable disk size(e.g. 
Mattson & Weaver 2004; Markowitz, Edelson & Vaughan, 2003). These arguments are based 
on temporal signatures; specifically the lack of a prompt response of the "reflected" emis- 
sion to the continuum assuming light travel times within reasonable accretion disk spatial 
scales. While an accretion disk torus has been suggested as an alternative, i.e. more remote, 
reflector, the same temporal signatures could be reproduced within the context of a central 
source and an ambient outflow. 
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We note the concerns of some authors, e.g. Miller et al. (2004), that the hot disk 
inner-disk component and high frequency QPOs ( few times 100 Hz) seen in accreting Black 
Hole (BH) at the highest inferred accretion rates would not be visible through an "optically 
thick" outflow. Those authors thus rule out any possibility of ambient spectral reprocessing 
and distortion of the iron line feature, such as broadening and red-wing enhancement, by 
the outflow. On the other hand, Laurent & Titarchuk (2004) infer the outflow Thomson 
optical depth r of the out-flowing medium from parameter fitting using XMM and ASCA 
measurements of broad, red-shifted iron lines. They find that r never exceeds 2-3 in 
any of the cases analyzed. Presumably, the observer sees the radiation of the BH central 
source through the "haze" of the moderate optical depth. Furthermore Titarchuk, Cui & 
Wood (2002), hereafter TCW02 give a precise model for the loss of the modulation due 
photon scattering. It follows from TCW02 that the outflow optical depth tt,o would need 
to be around 16 and higher in order to suppress QPO amplitude of frequency 100 Hz [see 
formula (5) in TCW02]. We have shown that this very optically thick outflow is ruled out by 
observations. Our results confirm the LT04 results in the sense that all effects of reprocessing 
and line distortion can plausibly occur in outflows characterized by moderate optical depth. 

Furthermore, recently Laming & Titarchuk (2004), hereafter LaT04, formulate a generic 
problem of the outflow illumination by the hard radiation of the central object and they 
calculate the outflow temperature and ionization balance as a function of the ionization pa- 
rameter. Natural assumptions regarding the X-ray spectral distribution of the central source 
radiation (Comptonization like spectra) and velocity distribution in the outflow (constant 
velocity wind) are applied. They find that iron K a photons are generated by the absorption 
of X-ray photons at energies higher than the K-edge (i.e > 7 keV) (see these K-edge features 
in Fig. 3). Electron scattering of the K a photons within the highly ionized expanding flow 
leads to a decrease of their energy (redshift) which is of first order in v/c (this is clearly 
illustrated in Fig.l) This photon redshift is an intrinsic property of any outflow for which 
divergence is positive. LaT04 find the range of the parameter Inner Radius/L40 (which is 
proportional to the inverse of the so called "ionization parameter" used in the literature) is 
about 10 13 cm and the range of the wind temperature is about 10 6 K when the observed K a 
lines are produced in the wind (where L40 is the source luminosity in 10 40 ergs s _1 ). They 
also find that the equivalent widths of red skewed Fe Ka originated in the wind is order of 
keV. 
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5. CONCLUSIONS 

We have developed an analytic formulation for the emergent spectrum resulting from 
photon diffusion in a spherically expanding Comptonizing media, characterized by two model 
parameters: an average number of scatterings in the medium N av and efficiency the energy 
loss in the divergent flow e. In this formulation of the Radiative Transfer problem, the number 
of photons emitted in the wind shell equals that which escape to the Earth observer (i.e. the 
photon number is conserved). Consequently, if the high energy photons lose their energy on 
the way out but the number of photons is conserved, it has to result in the accumulation of 
the photons at a particular lower energy band. 

Application of the model to high-energy spectra of several compact binaries, which 
have been previously modeled within the Compton reflection scenarios, seems to lead to a 
satisfactory representation of the data. We suggest that in some instances, the apparent 
~ 10— keV continuum enhancement seen in Galactic and extra- galactic sources is due to 
downscattering effects associated with an outflow rather than to reflection by a disk. 

We may also conclude that scattering and absorption of the primary line photons in 
the relatively "cold" outflow (that temperature is few times 10 6 K) lead to the downscattering 
modification of the continuum and to the formation of red-skewed lines (that is a more natural 
and probable mechanism than the general relativistic effects in the innermost part of the 
accretion flow). 

In future work, we will employ detailed Monte-Carlo calculations to further explore the 
range of validity. In addition to Galactic binaries, extra-galactic sources may be characterized 
by out-flowing plasma. With the expanding database of INTEGRAL observations, which 
will hopefully lead to an increasingly accurate characterization of the high-energy continua, 
it may be possible to to explore this idea. In addition to the continuum, we will explore the 
possibility of the line feature formation in the wind. 
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Table 1. Source Sample and Comptonization Model Parameters 



Source ID 


Observatory 


r 


N 


E* 


e 


Flux 




DoF 


Cyg X-l 


INTEGRAL 


1.74 


0.8 


175 


0.12 


10.01 


1.544 


1359 


Cyg X-l 


RXTE 


1.51 


1.7 


219 


0.12 


17.1 


1.651 


80 


GX339-4 


RXTE 


1.63 


3.2 


299 


0.10 


4.27 


0.988 


123 


GS1353-64 


RXTE 


1.35 


1.9 


103 


0.10 


5.2 


1.015 


122 


Cyg X-3 


INTEGRAL 


2.61 


1.6 


160 


0.10 


5.57 


1.218 


1550 


Cyg X-3 


RXTE 


2.03 


1.01 


270 


0.10 


13.1 


1.068 


68 



Flux is for 3-100 keV, units of 10 9 ergs cm 2 s 1 
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Fig. 1. — On the left side: Schematic diagram depicting the wind geometry. The outflow 
(wind) originates at the inner radius. The electron optical depth of the wind is of order 
unity. A photon emitted near the inner boundary and subsequently scattered by an electron 
moving with velocity v 1; impinges on an electron moving with velocity v 2 as shown. The 
change in frequency is f 2 = v\ [1 + (vi — v 2 ] • n/c) where n is a unit vector along the path of 
the photon at the scattering point. In a diverging flow (vi — v 2 ) • n/c < and photons are 
successively redshifted, until scattered to an observer at infinity. The color of photon path 
indicates the frequency shift in the rest frame of the receiver (electron or the Earth observer) . 
On the right side: In a converging flow (vi — v 2 ) ■ n/c > and photons are blueshifted. 
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Fig. 2. — a) ETe diagram representing the thermal Comptonization spectra that result 
from photon propagation through a divergent outflow. The photon source is in the center 
of the outflow shell. Solid line: the emission of the central source for which a spectrum is 
E~ a exp(— E/E*) where a = 0.5 and E* = 50 keV. Dashed, dot-dashed and dotted lines: the 
emission not escaping radiation from the outflow shell for that average number of scattering 
N av = 2, 3, 5 respectively. The effective coefficient of the outflow divergence e = 0.2. b) 
Similar to Figure 2a, except that a ratio of the Comptonization models to the exponentially 
cutoff power-law spectrum is plotted. 





Fig. 3. — a) Comparison between our Comptonizaton model (solid curve), and a simple cutoff 
power-law form, i.e. ~ KE~ r exp(—E/Eo) (dashed curve). Both curves resulted from fits to 
the overlayed data points, from GX 339-4. No additional parameters (such as absorption or 
emission due to Fe) were included. It is evident that the Comptonization model represents 
the data more accurately above ~ 10 keV (the improvement in xl is about 10%). As detailed 
in the text, we interpret the "excess" flux in that regime as being due to downscattering of 
hard photons in a diverging outflow, b) Fit residuals (for GX 339-4) above 10 keV for the 
cut-off powerlaw model (triangles) and our out-flow Comptonization model (filled circles). 
As evidenced from this figure, the downscattering effects provide a viable explanation of the 
10-100 keV "excess" continuum. 
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Fig. 4. — Example of the the application of our model to observational data. In this case, 
flaring state RXTE data for GS 1351-64 are shown. This source has been noted by others 
as a strong "Compton reflection" source among the Galactic X-ray Binaries. The count rate 
data,the folded model, and the residual are plotted. In this case, the chi-square per degree 
of freedom was xl — 1-02. 
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Fig. 5. — Similar to Figure 4, but the data and model fit are plotted in photon space. 



